Opto-Mechanical Pattern Formation in Cold Atoms 
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Transverse pattern formation in an optical cavity containing a cloud of cold two-level atoms is 
discussed. We show that density modulation becomes the dominant mechanism as the atomic tem- 
perature is reduced. Indeed, for low but achievable temperatures the internal degrees of freedom of 
the atoms can be neglected, and the system is well described by treating them as mobile dielectric 
particles. A linear stability analysis predicts the instability threshold and the spatial scale of the 
emergent pattern. Numerical simulations in one and two transverse dimensions confirm the insta- 
bility and predict honeycomb and hexagonal density structures, respectively, for the blue and red 
detuned cases. 
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Pattern formation can be defined as the spontaneous 
emergence of spatio-temporal structures in nonlinear 
systems driven far from equilibrium The spatial 

structure emerges from an initially homogeneous state 
as a consequence of the interplay between local non- 
linearity and spatial mechanisms such as diffusion and 
diffraction. In optical systems, spatial structures of the 
fight intensity are generated in the plane transverse to 
the direction of propagation after interaction with a 
nonlinear medium in the presence of feedback. Several 
geometries have been shown to provide the required 
feedback: two counterpropagating beams ring or 
Fabry-Perot cavities [3|-|5[, or the single mirror feedback 
arrangement [B-Q- Moreover, theoretical and experi- 
mental investigation demonstrated that different media 
can support the formation of transverse structures, 
such as (hot) atomic gases 0, H, 0, [3, Liquid Crystal 



11 



In 



Light Valves 10] or photorefractive media 
particular, atomic media present the appealing feature 
of being describable from first principles combining the 
microscopic equations for the medium and the Maxwell 
equations for the incident and generated radiation. It 
is well known that diffraction and optical nonlinearities 
can induce spatio-temporal structures inside the medium 
through modulation of populations and coherences. 
An additional mechanism for spontaneous self- 
organization appears in atomic media cold enough 
for optical forces to compete with thermal effects 
sufficiently strongly to modify the density distribution of 
the sample. Nonlinear optical effects involving the me- 
chanical effect of light on cold atoms are well known and 
have been predicted and demonstrated in cold atomic 
gases 12 -IJ] and Bose-Einstein Condensates 15, [l6| . 
In particular, the possibility of using such effects for 
spontaneous filamention has been proposed in [l7|. The 
nonlinear interplay between these opto-mechanical forces 
and the induced density modulations that they create 
and respond to has been shown to give rise to collective 
or cooperative light scattering 12h161 Il8l421| in atoms 
and self- focusing 122.], and four- wave mixing [2^ in 



solutions of dielectric spheres. However, the formation 
of opto-mechanical structures, arising from spontaneous 
symmetry breaking in the plane orthogonal to the 
pump axis due to opto-mechanical forces, has not been 
addressed in these studies, and is the main topic of this 
Letter. For counterpropagating beams, a lowering in the 
threshold for transverse self-organisation on the focusing 
side 
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of the nonlinearity has been proposed 
Experimental evidence of the formation of transverse 
structures in cold atomic media has also been found 
in 2^ 27 1 . These previous studies typically emphasize 



the interplay between mechanical density redistribution 
effects and the nonlinearities arising from the internal 
degrees of freedom of the atom, and often involve multi- 
ple optical beams forming wavelength-scale lattices, and 
perhaps also optical polarization effects related to the 
multi-level quantum structure of the atoms. 
In contrast, we consider a very simple, and hence general, 
system, of ground-state atoms interacting with a single 
coherent optical field. The linear dielectric response 
of the atoms, which is responsible for the refractive 
index of the cloud, means that the atoms will move up 
or down any transverse gradient in the optical field. 
In turn, the refractive effects of non-uniform atomic 
density will lead to phase gradients, and thus eventually 
intensity non-uniformities, in the optical field. As we will 
show, this simple mechanism readily produces positive 
feedback, and thus transverse instability, independently 
of the the sign of the atomic response. We illustrate 
this very general instability mechanism for the case of a 
Doppler cooled two-level atom cloud within a planar ring 
cavity driven by a monochromatic plane wave optical 
field. This simple fundamental atom-field coupling is, 
of course, already present in all the above-mentioned 
experiments, and in many of the models. Its impor- 
tance has not previously been apparent because of the 
complexity of the systems considered. We believe that 
the results of our simple model will be important for 
the interpretation of these more complex experiments 
and models, and hence for future progress in cold- atom 
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FIG. 1: (color online) A sample of two level atoms with 
thickness I and density A^o is inserted in a planar ring cavity 
of effective length L and laser-cooled at a temperature T. L 
can be controlled by adjustment of the intra-cavity lenses, 
and is therefore distinct from the pysical cavity length. A 
plane wave of amplitude y and frequency uj, detuned from the 
resonance ujat by A, pumps the medium in the cavity. The 
transmittivity of the mirrors is r. A typical configuration can 
be obtained exploiting the D2 line (A = 780 nm) in a sample of 
*^Rb of thickness/ = 10 mm and density Ao = 7x 10^°at/cm^, 
cooled to twice the Doppler temperature {Td = 150 ^K) and 
inserted in a cavity with mirror transmittivity r — 0.1. 



optics. 

We consider a sample of N identical, non interact- 
ing two-level atoms inserted in a planar ring cavity and 
cooled to temperature T (see Fig. [T|). We allow for 
a transverse redistribution of the sample density under 
opto-mechanical forces in the form N{x,t) = NQn{x,t), 
where A'o is the average density of the sample and 7i(x, t) 
encodes a spatial density modulation. Such a modula- 
tion will enter the expression for the sample suscepti- 
bility, with regions of higher density corresponding to 
larger responses of the medium to the incoming radi- 
ation. The material susceptibility is cast in the form 
X = N{x,t)xc{^,t), where N accounts for the den- 
sity redistribution effects and Xc represents the elec- 
tronic susceptibility of the medium. Adiabatic elimina- 
tion of populations and coherences in the optical Bloch 
equations of a two-level system are known to produce 
an intensity-dependent electronic response of the form 
Xe = Xo(|/(x, where / represents the amplitude of 
the electric field [4|. The medium susceptibility x then 
acts as a nonlinear source term in the field wave equa- 
tion, which can be written in the slowly varying envelope, 
rotating wave, paraxial and mean field approximations 
as i: 



f = -il + i9)f + y-^n- 
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Eq. ([T]) contains only adimensional quantities, with the 
field / rescaled to the saturation intensity at resonance 
and the time normalized to the cavity losses (see the 
term — /). A is the light-atom detuning (in units of 
the coherence decay rate r/2), while V5_ = 
denotes the transverse Laplacian and describes diffrac 
tion. The cavity imposes a linear phase shift governed 
by the cavity detuning 9, and is pumped by a plane 
wave of (normalized) amplitude y. Definitions are cho- 



92 



sen so that A > (A < 0) corresponds to blue (red) de- 
tuned beams and thus to self-focusing (defocusing) non- 
linearities. The strength of the complex susceptibility is 
7 = 2C(l-|-iA)/(l-|- A^), with absorption and dispersion 
captured by its real and imaginary parts, respectively. 
The cooperativity parameter C contains the dependence 
on the sample density A^o smd the mirror transmittivity 
T. Spatial coordinates are normalized to the diffraction 
length y/a = ^ \L/Atit, where A is the radiation wave- 
length and L the effective cavity length (see Fig. [1]) . Note 
that terms in the Maxwell equations varying as fi, h are 
neglected in deriving Eq. ([T]). 

If n is considered uniform, Eq. ([T]) becomes a standard 
model of cavity nonlinear optics, showing optical bista- 
bility and pattern formation for appropriate parameter 
choices. If the nonlinear term in Eq. ([T]) is neglected, no 
instability is possible for constant n, and the role of the 
atoms is simply to modify the cavity losses and resonance 
frequencies through the complex linear susceptibility. If 
we allow the density to enter as a dynamical variable, 
however, the system behavior becomes qualitatively dif- 
ferent. To proceed, it is necessary to specify an equation 
for the dynamics of the density modulation n, which is 
coupled back to the field through the action of optical 
forces. In the limit of large detuning scattering forces are 
negligible, and the sample is subject only to a conserva- 
tive dipole potential [/dip = (?irA|/|2)/4(l + A^ + |/|2). 
Assuming a strong viscous damping of the momentum 
distribution, e.g., due to the presence of an optical mo- 
lasses, a Fokker-Planck equation can be derived for the 
dynamics of the atomic density [111, 111] : 
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We recognise the first term on the right side of Eq. ([2]) 
as the divergence of a drift current originating from 
the transverse dipole forces, potentially leading to non- 
uniform density, while the second term, diffusion, op- 
poses such non-uniformities. The parameter a charac- 
terizes the relative strength of these opposing tenden- 
cies. Clearly density modulation is favored by large de- 
tuning (A) and, importantly, by low temperature. In 
the limit of high temperatures (cr — 0) diffusion drives 
the atomic distribution towards the homogeneous state 
(n = 1), and the model reproduces the results corre- 
sponding to hot atomic vapours. The stationary state 
for the density modulation is given by the equilibrium 
distribution [iB,!!!: 
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yexp(-[/dip/fcBr) 
Jy dx exp {-Udip/kBT) 



(3) 



Eqs. ([T]) and 1^ describe the coupled dynamics of the 
intra-cavity field and the two-level sample when both 
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electronic and opto- mechanical effects are present. This 
kind of system, and the role played by the temperature, 
has been addressed for the arrangement of two counter- 
propagating beams in 24 , [2^ . Here we study the much 
simpler situation of a unidirectional beam and, more- 
over, negligible electronic nonlinearity. We show that 
opto-mechanical effects alone are capable of providing a 
pattern-forming instability, through density redistribu- 
tion. For large detuning (|A| ^ 1) scattering forces and 
absorption are negligible compared to dipole forces and 
dispersion, respectively. If we also neglect the electronic 
nonlinearity, the system (IT][2]) reduces to 



f ^ -{1 + i9)f + y - i^n f + iVlf 
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where 7 — = 2CA/(1 + A^) accounts for lin- 

ear dispersion and nonlinear terms in -|- A^) have 

been neglected. We remark that this limit can be exper- 
imentally feasible. For a 10 mm thick sample of ^^Rb 
with density No = 7 x 10^" at/cm'^ at a temperature 
T = 300 /xK, interacting with a laser beam detuned by 
|A| ~ 100 linewidths from the D2 line and mirror trans- 
mittivity of 10%, for instance, one obtains C ~ 225 and 
\a\ ~ 25, and we will find that the electronic nonlin- 
earity is indeed small at the threshold for density-driven 
pattern formation. To demonstrate this, we perform a 
linear stability analysis of the system The flat, sta- 
tionary state of dH) is given by fs = y[l + i{9 + j)]^ , 
Ug = 1. We perturb this fiat solution as / = fs + 6f{x, t), 
n = 1 -|-(5n(x, t) and linearize the system Q. When look- 
ing for static instabilities we write the perturbations in 
the form ^ e^'^'^^e'^'^'^''* (with A real) and impose the con- 
dition of marginal stability, i.e. A(q) = 0. A threshold 
condition is found for the control parameter / = |/sp: 



9c,OM = l-{0 + l), 
^th — • 
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(5a) 
(5b) 



Here qc,oM denotes the critical wavenumber of the sys- 
tem (HI, i.e. the first transverse mode to become ex- 
cited when increasing the control parameter. The critical 
wavenumber can be controlled by varying the cavity de- 
tuning 6: we set 9 + J = —1 {9 ~ —5.5), which gives 
9c OM ~ 2- -^th^ represents the minimum value of inten- 
sity such that the growth rate crosses zero, A(qc,oM) — 0. 
We find ~ 88.9 for our choice of parameters. The 
low-excitation assumption is therefore confirmed since 
/t^^/(l + A2) ~ 0.009 < 1. 

Figure [2] shows the threshold curves I{q^) for our choice 
of parameters. Fig. [2^ compares the threshold for the 
purely opto-mechanical system of Eqs. (jlj (full line) with 
that of the full system involving both electronic and me- 
chanical effects, Eqs. (PH)), for blue (blue dashed line) 




FIG. 2: (color online) Left: pattern-formation threshold 
curves for Eqs. (|4|) (full line) and for the full system Eqs. (|ll2p 
for blue (dashed line) and red (dash-dotted line) detuning. 
|A| = 100, |cr| = 25, C = 225 and £» 7 = -1. Right: on a 
much larger intensity scale, thresholds for the purely mechan- 
ical system (full line), the complete system (blue dots) and 
the blue-detuned case of the saturable Kerr model (dashed 
line) without density modulation. 



and red (red dot-dashed line) detuning. Electronic non- 
linear effects clearly cause only small corrections to the 
threshold ([5]) from the purely opto-mechanical model. 
The minimum thresholds for the full model are found 
to be /th — 90.8 (blue detuning) and /th — 100.7 (red 
detuning). The critical wavenumbers are ~ 2.07 (blue 
detuning) and q'^ ~ 2.01 (red detuning). The instability 
behavior of the full system is clearly dominated by the 
opto-mechanical modulation of the atomic density, with 
electronic effects negligible in first approximation. Fur- 
ther evidence is displayed in Fig.[2]3, which compares the 
thresholds for the purely mechanical system of Eqs. ^ 
(full line) with that for a saturable Kerr medium with no 
density redistribution effects, i.e. a hot two-level medium 
(dashed line), which is more than two orders of magni- 
tude greater (for blue detuning, A — 100: there is no 
hot-atom instability for red detuning for our choice of 9). 

We note that instability thresholds arising from opto- 
mechanical effects are independent of the sign of the de- 
tuning, as a change in the sign of A (and thus 7) is com- 
pensated by a change in the sign of a. This is strongly 
reminiscent of the artificial Kerr media studied by Ashkin 
and collaborators 22| , where dielectric particles are sub- 
ject to light forces. A change in the sign of the parti- 
cle polarisability (i.e. from blue-detuned to red-detuned 
beams) is compensated by the fact that positively (neg- 
atively) polarised particles are pushed to the maxima 
(minima) of the field intensity. The value of the tem- 
perature is crucial in determining the interplay between 
the opto-mechanical and the electronic nonlinear mech- 
anisms. The parameter a represents in fact the ratio 
between the dipole energy h{uj — cjat) and the thermal 
energy ksT. By increasing this ratio, opto-mechanical 
nonlinear effects arising from dipole forces become dom- 
inant over the electronic nonlinearities from saturation 
of the involved transition. Internal degrees of freedom 
thus become negligible close to the instability threshold 
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FIG. 3: (color online) Intensity (black solid lines) and density 
(dashed lines) transverse structures in one transverse dimen- 
sion obtained from numerical simulations of the system ([iJS]) 
at about 3% above threshold. Parameters are as in Fig. [S] 
with the blue-detuned case on the left panel and the red de- 
tuned case on the right panel. Electronic nonlinear effects can 
be neglected, so that the system undergoes a pattern-forming 
bifurcation with critical wavenumber ~ 2.07 {q^ — 2.01) 
and threshold intensity Jth = 90.8 (Jth ~ 100.7) in the blue 
(red) detuned case. Density structures form according to the 
values of detuning and temperature, see Eq. ([3} ■ 




FIG. 4: Opto-mechanical hexagon formation at about 7% 
above threshold. Parameters are as in Fig. [21 with the blue- 
detuned case in the top panel and the red detuned case in 
the bottom panel. Field intensity is depicted on the left, the 
corresponding density modulation on the right. 



which can thus be interpreted as due to opto-mechanical 
effects alone (see Fig. [2]). We remark that, although \a\ 
can, in principle, always be increased to very high val- 
ues by lowering the temperature, our parameters do not 
require sub-Doppler temperatures, because our large de- 
tuning parameter weakens the electronic reponse relative 
to the opto-mechanical. 

Numerical simulations in one and two transverse di- 
mensions have been performed to test the predictions 
obtained from the linear analysis. The models ([l][2]) and 
(HI) have been integrated using a second-order Crank- 
Nicholson method, with the density dynamics given by 
^ . Periodic boundary conditions are imposed over a do- 
main of 5 critical wavelengths Ac = in/qc. The trans- 
verse domain is discretized using 112 points for one- 
dimensional (ID) simulations, and a square grid of 80 x 80 
points for two-dimensional (2D) simulations. Time is dis- 
cretized with step = 5 X 10""^. Fig. |3] presents results of 
ID numerical integration of the full system (IT][2]) for the 
parameters of Fig.[5]and intensities slightly above thresh- 
old, / = 92.6 (blue detuning) and / = 103.7 (red detun- 
ing). In Fig. [3] the dashed lines correspond to the density 
distribution: from Eq. ([3|) the absolute value of a governs 
the amplitude of the modulation, while its sign deter- 
mines whether atoms bunch in regions of high (cr < 0) or 
low (cr > 0) intensity, see Fig. |3l Fig. |4] presents results 
of numerical integration in 2D of the simple system ^ , 
for the same parameters as Fig. [21 for both signs of cr and 
I = 95, about 7% above threshold. The opto-mechanical 
nonlinear mechanism leads in both cases to the forma- 
tion of hexagonal structures, as is usual in systems with 
intensity-dependent cubic nonlinearities [1^. Note that 
the intensity pattern is bright hexagons in both cases, 
but the density structure is honeycomb- type for blue de- 
tuning. This is because the linear refractive index of the 
atom cloud is less than unity, so that a "hole" in the 



cloud has relatively high index, and so can guide light. 
We thus interpret the upper panels of Fig. [H as a self- 
organized hexagonal network of waveguides formed by 
the expulsion of atoms from the guided beams. Since the 
overlap of light and atoms is reduced in the fully-formed 
pattern the opto-mechanical nonlinearity has a strong ef- 
fective saturation, and this pattern is very stable. For red 
detuning, in contrast, atoms and light attract each other, 
so that the light is now guided by atomic filaments of high 
density (lower panels Fig. [31). As evidenced by the varia- 
tion of the amplitude of the peaks in Fig. 31 (lower panel), 
here the interaction is enhanced by the pattern forma- 
tion, and it is perhaps unsurprising that the stability of 
the pattern is much poorer for red detuning, though de- 
tailed investigation of stability issues is beyond the scope 
of this Letter. 

With a view to experimental observation of opto- 
mechanical pattern-forming instabilities, we note that an 
effective cavity length L of a few centimeters would lead 
to pattern scales of the order of Ac ^ 100 //m, requir- 
ing beam diameters of at least ~ 500 /im for the de- 
velopment of a well-defined transverse structure. For 
standard Rb molasses at temperatures of T ^ 100 /xK 
the threshold intensity is of order lOOmW/cm-^, so that 
2 mW of intra-cavity power should be ample to generate 
opto-mechanical hexagons. 

We have investigated a mechanism for transverse op- 
tical pattern formation dominated by density modula- 
tion effects due to dipole forces exerted by light on a 
cold atomic medium. By considering a very simple ring- 
cavity configuration we have been able to isolate and 
highlight the role and importance of opto-mechanical ef- 
fects in light-atom interaction at low temperatures. The 
threshold and the spatial scales of the resulting patterns 
have been confirmed by numerical simulations in one and 
two transverse dimensions. The required atomic temper- 
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aturcs and optical intensities are well within experimen- 
tal capabilities, and indeed we suggest that transverse 
density modulation effects may well be present in previ- 
ous studies involving more complex configiirations. Our 
simulations predict the formation of self-organized hexag- 
onal wave-guiding filaments, the filaments being atom- 
poor or atom-rich depending on the optical frequency. 
The resulting coupled light-matter structure can be inter- 
preted as a self-organized and self-loaded optical lattice, 
optical lattices being the workhorse for applications of 
cold atoms solid-state physics and quantum information. 
As for other systems which form hexagonal patterns, we 
would expect the existence of stable single filaments, i.e. 
localized states or dissipative solitons. 
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